Wakefield generation in magnetized plasmas 
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We consider wakefield generation in plasmas by electromagnetic pulses propagating perpendicular to a strong 
magnetic field, in the regime where the electron cyclotron frequency is equal to or larger than the plasma 
frequency. PIC-simulations reveal that for moderate magnetic field strengths previous results are re-produced, 
and the wakefield wavenumber spectrum has a clear peak at the inverse skin depth. However, when the cyclotron 
frequency is significantly larger than the plasma frequency, the wakefield spectrum becomes broad-band, and 
simultaneously the loss rate of the driving pulse is much enhanced. A set of equations for the scalar and vector 
potentials reproducing these results are derived, using only the assumption of a weakly nonlinear interaction. 

PACS numbers: PACS: 52.35.Mw, 52.40.Db, 52.40.Nk 



I. INTRODUCTION 

Since the pioneering work in Ref. [IB, where the excita- 
tion of a plasma oscillations by a short laser pulse was con- 
sidered, much interest has been devoted to wakefield gen- 
eration in plasmas. To a large extent this has been due to 
the wakefield acceleration concept (see, e.g., 01-0]), where 
the longitudinal field traveling close to the speed of light in 
vacuum is used as an effective source for particle accelera- 
tion. However, the wakefield properties has also been study 
from a more theoretical point of view 1I7I-H2I1. Most studies 
of wakefield generation has been done for unmagnetized plas- 
mas. This is justified by the fact that in most experiments 
(O c j <X)p <C 1, where <o c = qBo/m is the electron cyclotron fre- 
quency, (O p = (noq 2 1 Cora) 1 / 2 is the plasma frequency Bo is 
the unperturbed magnetic field strength, uq is the unperturbed 
electron density, and q and m is the electron charge and mass 
respectively. As can be expected, the wakefield properties is 
not much affected by the magnetic field in this regime. How- 
ever, for C0c/<X)p > 1, and for the exciting source propagating 
non-parallel to the the external magnetic field, the wakefield 
is significantly altered by the non-zero Bo 1 10-13:!]. In partic- 
ular, as found by e.g. Refs. lfl2l U3I1 . for the EM-pulse trav- 
eling perpendicular to the magnetic field, the excited mode 
becomes an extra-ordinary mode, with a significant electro- 
magnetic part, and a non-zero group velocity. 

In the present paper we will re-consider the excitation of 
the extra-ordinary mode wakefield in a strongly magnetized 
plasma, allowing for (O c /(O p 3> 1. In this regime the group- 
velocity of the wakefield may approach the speed of light in 
vacuum c. In then turns out that the excitation process is sig- 
nificantly altered, as compared to the regime of a more modest 
value of Bo- From Particle-In-Cell (PIC) simulations we find 
that when Bo is increased, the wakefield changes from an ap- 
proximately monochromatic field at low (o c / (O p , to a broad- 
band one at (Q c /(O p S> 1. At the same time the energy loss 
rate of the exciting high-frequency ordinary mode is much in- 
creased with increasing co c /co p . A reduced system of equa- 
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tions is derived, generalizing previous Eqs. lfl2ll avoiding the 
division into fast and slow time scales. The reduced equations 
are solved numerically, and a good agreement with the full ID 
PIC simulations is found. 



n. PIC-SIMULATIONS OF WAKE FIELD GENERATION 
IN A MAGNETIZED PLASMA 

Here we will study an ordinary mode propagating parallel 
to an external magnetic field, Bo = Box. In the absence of the 
magnetic field, it is well known that the ponderomotive force 
due to a short electromagnetic pulse will excite a wakefield of 
electrostatic oscillations lll|2t]. Here we will focus on the ef- 
fects due to the external magnetic field lfl2l [l3ll . We note that 
this is more likelyto be relevant for experiments in the mi- 
crowave regime [14], rather than the regime of optical lasers. 

The ID Particle-In-Cell simulation (LPIC++)||l5t] is carried 
out to study the effect of external magnetic field on the wake- 
field generation in magnetized plasma. The typical simulation 
geometry is shown in Fig. 1, where the exciting electromag- 
netic (EM) high-frequency pulse is taken to be an ordinary 
mode. For all the results presented here the EM field ampli- 
tude of the ordinary mode driver is considered as 0.2 (in di- 
mensionless unit «o = eE/mcoc, with CO being the frequency), 
space and time are taken in units of laser wavelength (A) and 
one laser cycle z = A/c respectively. The 800 nm laser with 
FWHM duration of 8t and an initially Gaussian shape prop- 
agates through the plasma along the z-direction with electric 
field along x and the laser magnetic field along the y-direction 
lflr3l . We have modified the code in order to extend its ability 
to incorporate the presence of an external magnetic field. The 
external magnetic field is taken to be along the jc-direction. 
The plasma of length 800 A is considered with a 40A of the 
ramp region at the front. In the ramp region the plasma den- 
sity increases linearly from to 0.000625 n c (Fig. [TJ, where 
n c = S{)mco 2 /q 2 . 

Due to the external magnetic field the wakefield gets both a 
longitudinal and a transverse part. In Fig. |2]the longitudinal 
part of the electric field wakefield is shown when varying the 
value of the external magnetic field. Similarly the transverse 
component is shown in Fig. [3] 

A number of features from the PIC-simulations agree with 
the model and results presented in Ref. II 1211 (see also the dis- 
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FIG. 1. Simulation geometry for the study. 
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FIG. 2. Spatial profiles of the longitudinal wakefield at 600t in pres- 
ence of external magnetic field with ratio CO c /co p = 0.5 (a), 2 (b), 4 
(c) and 8 (d) respectively are presented using PIC simulation. 
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FIG. 3. Spatial profiles of the transverse wakefield at 600t in pres- 
ence of external magnetic field with ratio CO c /(Op = 0.5 (a), 2 (b), 4 
(c) and 8 (d) respectively are presented using PIC simulation. 



tirely accurate. Furthermore, the theory of Ref. [fPH be- 
comes increasingly inadequate when the magnetic field 
strength is increased such that (O c /(O p 3> 1. In partic- 
ular, Ref. [12] predicts that the loss-rate of the driv- 
ing pulse is approximately independent of G) c /o) p . In 
Fig. g]the total energy °c S™ F + \E%\, where Sf F = 
\E x B y — E x B y — (B x — Bo)E y \ is the EM part, contained 
in the wakefield is plotted as a function of time for dif- 
ferent values of co c / C0 p . We see that the loss rate of the 
driving pulse due to the generated wakefield is approx- 
imately constant until co c / CO p > 2, but for larger values 
of this ratio it increases rapidly. 



cussion and simulation results in Ref. lfl3lo 

1. The ratio of the electromagnetic to electrostatic field 
field amplitude scales as co c / CO p when Bo is varied. This 
is in accordance with the wake field being an extra- 
ordinary wave mode. 

2. By observing the position of the drop in the wake field 
amplitude, we can deduce that wake field propagates 
with a group velocity v^ w f« ceo, 2 /((O p + 03^), which is 
in accordance with Refs. IU2l[l3n 

However, there are also some properties of the PIC- 
simulations that contain new results, as compared to previous 
findings. 

1. The wakefield does have a small but nonzero trailing 
part at a distance more than (v g — Vg W f) T, where v g « c 
is the group velocity of the ordinary mode, and T is the 
time after the entrance of the driving pulse. The relative 
energy content of the wakefield in this tail is larger for 
larger value of co c / co p . 

2. The approximation that the evolution of the wakefield 
is slow in a frame moving with the pulse is not en- 
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FIG. 4. Temporal evolution of the total energy of the wakefield (ES 
and EM components) for different values of the external magnetic 
fields using PIC simulation. 

One way to understand the deviance from the models in 
Ref. lfl2tl is to investigate the validity of the basic assump- 
tions. A key ingredient in the derivation is that the evolution 
should be slow in a frame moving with the group velocity 
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of the driver. A consequence of this is that the phase veloc- 
ity of the wakefield should match the group velocity of the 
driver. Since the wakefield is dispersive, such a condition de- 
termines the dominating frequency and wavenumber to be ex- 
cited. Noting that the plasma is underdense with respect to the 
driver frequency, the group velocity of the driver is close to c, 
which matches the wakefield phase velocity in case the ex- 
cited wakefield frequency is co p and wavenumber co p /c. Such 
a condition can be seen to hold for the leading part of the 
wakefield, but for a strong static field this part actually makes 
up for a small part of the energy content of the excited extra- 
ordinary mode. This is illustrated in Fig. where the fre- 
quency spectrum of the wakefield is computed as a function 
of the static magnetic field for a fixed position and time of 
interaction. For (Q C /(0 P > 2 we see that the spectrum ceases 
to peak at co p , and as a consequence the assumption that the 
evolution of the wakefield is slow in a system moving with 
the group velocity cannot be accurate. In the rest of the paper 
the aim is to investigate the physics of wakefield generation in 
the regime where the basic equations of Ref. [ 12] needs to be 
improved. 



III. THE WEAKLY RELATIVISTIC 1-D MODEL 

For ID-perturbations (with variations along z), the cold rel- 
ativistic fluid equations can be written. 



and 



dt +Vz dz 



dn d , . 



q \ q d<t> 
yv-— A =— — z 
m J m az 



C0 C \ x x 



(1) 



(2) 



where we have assumed the field to be of the form E = 
-dA/dt — V4>, B = V x A + Bo, the constant magnetic field 
is Bo = Bqx, 7= (1 — v 2 /c 2 )'/ 2 is the relativistic factor, and 
fflfc = qBo/m is the (constant) cyclotron frequency. For the 
case where the vector potential was a weakly modulated func- 
tion, A =A(z,r)exp(fe — (Dt) +c.c, where c.c. stands for 
complex conjugate, Ref. 1U2I1 derived a coupled set for A (with 
A=Ax) and 4>, where A was the vector potential amplitude for 
the ordinary mode of high-frequency (o» co p ), and <t> de- 
scribed the electrostatic potential of the extra-ordinary mode. 
It was also noted in Ref. lfl2ll that the wake-field was also 
found to have en electromagnetic part, that could be computed 
from <t>. The main modification of the wake-field generation 
due to the static magnetic field Bo was to introduce an elec- 
tromagnetic component of the wake field, which in turn re- 
sulted in a non-zero group velocity. The PIC-simulations of 
the previous section confirms that these properties are essen- 
tial ingredients when the external magnetic field is introduced. 
However, when the ratio co c / C0 p is increased beyond unity, we 
also saw that additional features were introduced. In particu- 
lar, the assumption of Ref. [12] that the fields evolve slowly in 
a frame moving with the group velocity, becomes less accurate 
with increasing value of co c /co p . In the rest of this section we 



will derive the coupled equations for an extra-ordinary mode 
driven by an ordinary mode, without any additional assump- 
tions than those of an ID geometry and a weak nonlinear- 
ity. In particular, we will avoid making a WKB-ansatz for the 
high-frequency ordinary mode. We use the Coulomb gauge 
and write the vector potential as A =Ax+A w y. Thus the com- 
ponent A w is associated with the electromagnetic component 
of the nonlinearly driven extra-ordinary mode, and <t> with its 
electrostatic component. The fields of the ordinary mode is 
derived solely from A. 

From the ^-component of the momentum equation (f2]i we 
obtain 



qA ( q 2 A 2 

Vr= 1 



2m 2 c 2 



which we substitute into Ampere's law, to give 
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2m 2 c 4 
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where co p = (n§q 2 j £§m) l l 2 is the plasma frequency associated 
with the unperturbed electron density no- Next we study the 
z-component of the momentum equation. We keep variables 
that are linear in the extra-ordinary wave-mode, and include 
quadratic nonlinearities in the ordinary wave mode variables. 
Essentially this nonlinearity turns out to be the ponderomo- 
tive force, but keeping also the second-harmonic term in ad- 
dition to the low-frequency contribution. The z-component 
of the velocity is then expressed in terms of the electrostatic 
potential (through the z-component of Ampere's law) and the 
ordinary variables are expressed in terms of A. Finally the y- 
component of the velocity is expressed in terms of A w (from 
the j-component of Ampere's law). The equation is then ex- 
pressed solely in terms of the potentials, and we obtain 



dt 2 



CO' 



<94> 
~dz~ 



c 



dt 2 



dz 2 



A w = - 



C0 2 q dA 2 
2mc 2 dz 



(5) 

Finally, we study the y-component of the momentum equa- 
tion, and substitute the velocities in terms of the potentials in 
the same way as previously, which leads to 



dt 2 



dz 2 



cot 



<94> 

Cfflby- (6) 



By using Poisson's eq. the density perturbation can be ex- 
pressed in terms of the potential, d 2 <P/dz 2 = qSn/eo, and 
Eqs. (lUl-© forms a closed set for A,A W and <t>. For zero 
external magnetic field co c — > 0, the coupling to the electro- 
magnetic part of the wake-field A w vanishes, and together 
with (0 form a closed set for A and <t>. By eliminating d<t>/dz 
in terms of A w using (O in Eqs.© and (0 we may get a cou- 
pled set of two equations for A and A w , which however has 
fourth order derivatives from the wave operator of the extra - 
ordinary mode. Specifically, we obtain 

'*L+rA 2 ^L 

. dt 2 + p \ dt 2 C dz 2 



co: 



-co; 



dt 2 



dz 2 



A w — 



C0 2 CO c q dA 2 
2mc 2 dz 



(7) 



4 



which confirms that the electromagnetic part A w becomes pro- 
portional to the external magnetic field through co c - 

In order to perform a numerical analysis, let us first put 
the weakly relativistic ID-model in dimensionless form. In- 
troducing the normalized variables ?„ = co p t, z„ = zCO p /c, 
A n = qA/mc, <t>„ = q<P/mc 2 , A mv = qA w /mc, we obtain 



and 



d 2 A 
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CO, d<t> 
A w = ^— (10) 
(Op dz 



where we have left out the index n denoting normalized vari- 
ables for notational convenience. It is straightforward to show 
that the above system conserves total wave energy [Wdz if 
the energy flux across the boundaries vanishes. The energy 
density W = + W w f where the wake field contribution to 
the energy density is 



W w{ = 



1 
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dzdt 



~dz~) 

2 
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dt 
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with the terms form left to right representing longitudinal 
electric field energy density, magnetic field energy density, 
perpendicular electric field energy density, longitudinal ki- 
netic energy density, and perpendicular kinetic energy density. 
Similarly, the energy density due to the driving ordinary mode 
is 



W d = 



1 



dA 
~dl 



dA 
Hi 



1 - 



d 2 <S> 

dz 2 



(12) 

where the terms from left to right represent electric field en- 
ergy density, magnetic field energy density, lowest order (lin- 
ear) kinetic energy density, density perturbation correction to 
kinetic energy density, and finally relativistic correction to ki- 
netic energy density. The energy conservation law for Eqs. 
(O- ( [Tol l together with the expressions (fTTT > and (fT2l will be 
used in the next section for calculating the energy loss rate of 
the driving pulse due to the wakefield generation. 



IV. NUMERICAL COMPARISON 

Next we will compare the weakly relativistic ID model 
with the results from the PIC-simulations. The temporal evo- 
lution of the quantities A(z,t), A w (z,t) and d<S>(z,t)/dz is 
studied by numerically solving the dimensionless Eqs. dS}- 
( fT0T > using a standard explicit Leap-Frog scheme combined 
with the following initial and boundary conditions, 

A w (z,0)=0, d<P(z,0)/dz = 0, A(z,0)=0 



A(0,f) = A cos(a)?)exp 



-4 



Tfwhm 



Here Tfwhm is FWHM pulse duration of the laser which is 
considered to be Tfwhm = 1-25 in dimensionless units and is 
equivalent to the laser pulse duration used in PIC simulations. 

The numerical solutions of Eqs. ([8]l- ( fTOb for the longitudi- 
nal parts of the electric fields are shown in Fig. |6]for differ- 
ent ratios of co c / C0 p . Comparing these curves with the results 
from the full simulations, shown in Fig. [2] we see that the 
agreement is excellent. Calculating the energy evolution of 
the wakefield (as based on the energy densities given in (fTTT l 
and (fT2l). which is displayed in Fig. [7] we furthermore get 
convincing agreement with the PIC-results shown in Fig. [4] 
A more direct quantitative comparison of the energy loss rate 
is also shown in Fig. [8] where the PIC-simulation and the 
numerical results from (f8~t-(fT0l are simultaneously displayed. 
We thus conclude that the wakefield generation can be seen as 
the extra-ordinary mode driven by the ponderomotive force, in 
agreement with the model (0- ( fTOb . However, in contrast to 
an assumption commonly made in connection with wakefield 
generation, the evolution is not necessarily slow in a frame 
moving with the group velocity of the driver. To be more spe- 
cific - for a low to modest ratio of co c / CO p , the assumption of 
slow evolution in the co-moving frame is adequate. However, 
starting roughly from the value 2 < to c / C0 P it is not appropri- 
ate to consider the evolution as slow, which is reflected by the 
broad-band nature of the wakefield spectrum (see Fig. |5). It 
should be noted that besides abandoning the slow evolution 
assumption when deriving the Eqs. dS)- (fTOb we have also re- 
frained from the WKB-approximation of the high-frequency 
driver, that were made in Ref. [12]. By contrast to the slow 
evolution assumption, the validity of the WKB-approximation 
for the driver does not depend only on the ratio of co c /co p . In- 
stead the parameter kL (where L is the pulse length, and k the 
high frequency wavenumber) naturally plays a major role. In 
Fig. [5] we see that the wakefield frequency spectrum extends 
almost a factor 6 above the value co p . For efficient wakefield 
generation typically co p L/c ~ 1. Thus, in case kL ;§> 6, the 
high-frequency part of the wakefield spectrum is still much 
lower than the frequency of the high-frequency driver, and a 
WKB-approximation for this variable would be valid, whereas 
for kL ~ 6 it should be avoided due to the interaction with 
the turbulent wakefield spectrum which has a component of 
comparatively high frequency. The parameters used in our 
case, kL — 40, lies somewhere in the intermediate range, and 
therefore we have taken the safe approach and avoided the 
WKB-approximation for the high-frequency driver, although 
it should clearly be applicable for a sufficiently large value of 
kL, even for the highest values of co c / C0 P that has been consid- 
ered in our paper. 



V. SUMMARY AND CONCLUSION 

In the present paper we have re-considered ID wakefield 
generation in magnetized plasmas. Based on ID PIC Simula- 
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FIG. 5. Fourier spectrum of the transverse wakefield calculated at 
595A in presence of external magnetic field with ratio (O c /o) p = 0.5 
(a), 2 (b), 4 (c) and 8 (d) respectively are presented using PIC simu- 
lation. 
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FIG. 7. Temporal evolution of the total energy of the wakefield (Eq. 
Jilt is integrated over space for each time step) for different val- 
ues of the external magnetic fields by numerically solving the set of 
equations. 
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FIG. 6. Spatial profiles of the longitudinal wakefield at 600t in pres- 
ence of external magnetic field with ratio CO c /co p = 0.5 (a), 2 (b), 4 
(c) and 8 (d) respectively are presented by numerically solving the 
Eqs. [U-Q0] Here we have converted the units of space and time in 
order to compare with the PIC results. 



tions we have deduced that previous results need to be mod- 
ified in the regime of large magnetic fields. For moderate 
magnetic field strengths previous results are confirmed, and 
the wakefield consists of an extra-ordinary wave mode with 
a dominating frequency (o p and wavenumber CO p /c. However, 
monitoring the dependence as a function of the frequency ratio 
(Oc/ CO p , we find that there is a gradual transition in the regime 
2 S @c/®p ;$ 4 where the wakefield spectrum of the extra- 
ordinary mode becomes broadband, and the loss rate of the 
high-frequency driver is enhanced.. 

Although we have deduced that wakefield generation is 
governed by Eqs. (0- ( [Tol l, we have not so far discussed why 
the assumption of a slow evolution in the co-moving frame 
must be abandoned for a large value of co c /co p . However, 
an explanation of this can be given in the form of a reductio 




FIG. 8. Energy loss rate as a function of ratio (O c /o) p is presented for 
PIC as well as analytical results. 



ad absurdum reasoning. Let us assume that evolution of the 
extra-ordinary mode is indeed slow in the co-moving frame. 
In that case, there must be a dominant wavenumber of the 
wakefield corresponding to the property v& w f — v g « c, where 
v g is the group velocity of the high-frequency driving pulse 
and VQ W f is the phase velocity of the wakefield. However, 
from the dispersion relation of the extra-ordinary mode we 
find that the group velocity for this wavenumber scales ac- 
cording to v g — > c for (O c / (Op — > °°. This behavior is to a large 
extent reflected in Figs|2][3]and|5] where the dominant mode 
follows the driver closely for the stronger values of the mag- 
netic field. On the other hand provided the wakefield should 
follow the driver closely (as it must when v g —> c) and then 
drop in amplitude suddenly, it cannot be well represented by 
a small spectrum surrounding a leading wave-number. Hence 
we must have broad-band wakefield spectrum for strong mag- 
netic fields, in which case the assumption of a slow evolution 
of the wakefield will cease to be adequate. 
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To some extent the break-down of the slow evolution of the 
wakefield is related to the initial conditions. After long inter- 
action times the system may gradually set up a wakefield with 
a dominating wave number, and at the same time the slow 
evolution regime can be approached. However, the stronger 
the magnetic field is, the longer time this will take. For a 
limited size of the plasma region and a strong magnetic field, 
in practice the slow evolution regime will never be reached. 
A practical consequence of this is that the loss rate of the 
high-frequency driver will be much enhanced, due to the ex- 
citation of the broad-band spectrum. Furthermore, the energy 
loss will be accompanied by a frequency decrease and a de- 
acceleration of the driver fl7||. 



So far only the basic consequences of a strong magnetic 
field has been explored. A more complete study, e.g. involv- 
ing the effects of varying the self-nonlinearity and the long- 
time evolution, is a project for further study. 
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